Experimental determination of nucleation scaling law for small charged particles 
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We investigated the nucleation process at the molecular level. Controlled sticking of individual 
atoms onto mass selected clusters over a wide mass range has been carried out for the first time. 
We measured the absolute unimolecular nucleation cross sections of cationic sodium clusters Na„ 
in the range n=25-200 at several collision energies. The widely used hard sphere approximation 
clearly fails for small sizes: Not only vapor-to-liquid nucleation theories should be modified but also, 
through the microreversibility principle, cluster decay rate statistical models. 

PACS numbers: 36.40.-c, 34.10.+X, 82.60.Nh 



The formation andgrowth of droplets, clusters and, 
more generally, very small particles is of considerable in- 
terest. Indeed these phenomena are essential for instance 
in clouds formation or crystal nucleation. The formation 
and growth of small particles can occur in different ways. 
The formation can for example start from a seed (het- 
erogenous nucleation), in the presence of a buffer gas or 
simply occur in a pure vapor (homogenous nucleation). 
The growth is, for the smallest sizes, essentially governed 
by successive attachment of single units whereas for big- 
ger sizes coalescence dominates. A huge amount of pa- 
pers have been devoted to homogeneous nucleation, with 
the aim of getting a correct quantitative theory of the 
phenomenon. This aim has not been reached yet. A 
large majority of these theories originates from the so- 
called "Classical Nucleation Theory" (CNT), introduced 
by Becker and Doring It was originally developed 
for neutral droplets but can be applied without funda- 
mental changes to ions (it is called in this case Classical 
Ion-Induced Nucleation Theory, or CUNT) [2, l3( . Ion- 
induced nucleation often predominates in nature because 
ions have a higher sticking cross section (SCS) than neu- 
tral particles. It has recently been demonstrated in the 
case of cirrus clouds condensation 0]. One crucial pa- 
rameter in nucleation is the growth rate, which depends 
on the size, defined as the probability for a particle to 
grow by one unit per unit time. It depends on the atomic 
flux and the SCS 0, [||. For large enough particles, the 
SCS is given to an excellent approximation by the hard 
sphere model. In the cluster size range however, i.e. for 
very small particles made of tenth to hundredth units, 
this approximation is likely to fail. Furthermore all nu- 
cleation theories based on CNT define a so-called critical 
size whose properties are crucial QH H, Hi . Unfortunately, 
in many syste ms ( water [H, [j| , metals jld| | , organic com- 



pounds [Tl, 13, H|) the critical size lies in the cluster size 
range where the hard sphere approximation may not be 
relevant. 

SCS's are also key parameters in the analysis of clus- 
ter decay rate. As a matter of fact, in statistical models 
based on the microreversibility principle (or detailed bal- 



ance), the decay rate is proportional to the SCS [14L Il5l] . 
Cluster dissociation energies as deduced from their evap- 
oration rates [l6| are thus influenced by sticking cross 
sections. 

The size dependence of the sticking cross-section has 
already been addressed theoretically (see e.g. 17, Gil). 



However, while many experiments have been devoted to 
the measurement of clusters evaporation rates (l6| . ho- 
mogeneous SCS's of mass selected clusters have never 
been measured to our knowledge. We present here the 
first extensive measurement of clusters SCS's as a func- 
tion of their size and collision energy. The experimen- 
tal setup has been described in details elsewhere [l9j ]. 
Sodium clusters are produced in a gas aggregation source. 
They are ionized by a hollow cathode discharge. They 
are then thermalized in a heat bath; their temperature 
Tj can be varied from 150 to 500 K. During cross sections 
measurement, Tj ~ 150 K in order to avoid evaporation 
at least until final products reach the detector. Wiley 
Mc-Laren type electrodes operate a first mass selection. 
Next, clusters enter an electrostatic device that reduces 
the kinetic energy spread of the mass selected clusters. 
It allows to decrease the kinetic energy to a few eV. 
Thermalized mass selected clusters then propagate slowly 
across a cell containing a sodium vapor (~ 10~ 4 mbar for 
T=200 °C). Finally, a second mass selection determines 
the new size distribution at the output of the cell. 

Our raw data are the relative number of clusters which 
undergo a sticking I/Iq. This quantity is measured as a 
function of the cluster size and kinetic energy. Assuming 
a Boltzmann velocity distribution for the atoms, the SCS 
ct„ is given by: 
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with a = Ef./ (nksT) and I is the length of the cell, p the 
atomic density inside the cell, n the size of the incoming 
cluster, T the temperature of the vapor and erf is the 
error function. Ek is the kinetic energy of the cluster 
in the laboratory frame (in the following, kinetic energy 
refers to the energy in the laboratory frame and colli- 
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sion energy to the energy in the Center-of-mass frame). 
The density of atoms in the cell p is deduced from the 
cell temperature through the vapor pressure curve. We 
carefully measured the temperature at different positions 
in the cell in order to determine as accurately as pos- 
sible the density of atoms from the lowest encountered 
temperature. The density of atoms in the cell remains 
nevertheless the main source of uncertainties. An error 
of 5 K on the cell temperature leads to a relative error 
of 25% on the determination of the cross section. These 
uncertainties can shift the curves collectively, but do not 
change the size dependence of the cross-section. 

Note that under our experimental conditions evapora- 
tion is completely negligible. We checked for evaporation 
effects on a few sample masses: For Na^ 10 we varied Ti 
from 150 to 250 K and we did not see any evaporation 
effect. This is not surprising since the lifetime for the 
product Na^ n is still 450 ms at 250 K. For Na^ we var- 
ied Ti from 150 to 200 K and found no change in the 
SCS: the lifetime for Na+ is still 63 ms at 200 K. For 
NaJ-g we varied Ti from 150 to 340 K and we did see a 
change in the SCS after a rather long plateau. Evapora- 
tion starts to affect the cross section at about 290 K. We 
find a lifetime of 300 /is for the product Na^ 6 at this tem- 
perature (Ec=20 eV) whereas the time of flight is about 
130 /is. It is for sizes smaller than n = 25 at Ek=20 eV 
and Ti=150 K that the lifetime becomes comparable to 
the time of flight. We can thus safely exclude the effect 
of evaporation for the size range presented in our results. 

Let us emphasize that contrary to previous experi- 
ments we work in the single collision regime with very 
well defined initial conditions (size, temperature and col- 
lision energy) and that we explored a wide range of sizes 
(n = 25 — 200). In the few studies where reaction cross 
section of mass selected particles are measured as a func- 
tion of size or collision energy reactants underwent evap- 
oration and/or the size range was quite small 2(J 21, 22| . 

SCS's have been measured at two kinetic energies. Fig- 
ure Q] presents our experimental results for Ek=10 eV and 
20 eV. The most noticeable feature is that the two curves 
are different for small sizes and merge at about n = 80. 
From this size on, they follow the geometrical scaling law 
in n 2 / 3 predicted by the hard sphere model (Figure [T]): 
a geo = 7ri? 2 n 2 / 3 , with R = 2.4 A (the Wigner-Seitz ra- 
dius for sodium is 2.1 A). For smaller sizes, the cross 
section departs from the hard sphere model. It drops be- 
low the geometric cross section for Ek = 20 eV whereas 
it increases for 10 eV. For sizes smaller than 80, the cross 
section depends on the collision kinetic energy: it de- 
creases for increasing kinetic energy. Figure [1] also dis- 
plays cross sections calculated in the frame of a Langevin 
model pel [23| . The cluster-atom interaction is described 
by charge-induced dipole potentials. Two potentials are 
involved, corresponding respectively to A r a+ + Na and 
Na n + Na + (see Fig. [5]) . These two states are asymp- 
totically separated by the ionization potential difference 
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FIG. 1: Log-Log plot of the SCS as a function of cluster size 
at two kinetic energies. Experimental results are compared to 
a Langevin model and a hard sphere model. The calculated 
mean cross sections in a vapor at 400 K are also shown. 



AIP between atom and cluster. In this frame, the stick- 
ing can be described as a harpooning mechanism: due to 
the much larger polarisability of the cluster compared to 
the atom, an electron transfer from the atom onto the 
cluster is energetically favored. The adiabatic potentials 
presented in Figure [2] are calculated using as a coupling 



the hopping integral in Na^ 24|. From the Landau- 
Zener formula, the propagation is found to take place on 
the ground state potential curve. In the Langevin model, 
the largest reactive impact parameter b max is reached 
when the maximum of the effective potential is equal 
to the collision energy. The cross-section is then given 
by On — fdm al . It is taken as geometrical when the 
Langevin cross section becomes smaller than the geo- 
metrical one. This simple model gives satisfactory re- 
sults for Ek = 10 eV. However, at Ek — 20 eV, cal- 
culated values disagree with the experiment: the model 
is unable to account for experimental cross sections be- 
ing clearly below the geometrical ones. More elaborated 
models exist [25], |2g, 127], |28|, |29j but they would also give 
a cross-section larger than the geometrical one. Moreover 



we have a charge- induced dipole interaction whereas 25 1 
deals with neutral particles and [26j consider molecules 
with permanent dipole. 

One possible explanation for this behavior is the pres- 
ence of non reactive excited states of the cluster. Indeed 
since the experiment is done at fixed kinetic energies in 
the laboratory frame, the collision energy E c increases as 
the size n decreases: E c = + 3fc f T -7-r. Non reactive 

c n+l 2 n+1 

excited states that can be reached as the collision energy 
increases are likely to be responsible for the drop in the 
cross sections for small n. Preliminary calculations based 
on a two-center jellium model show the appearance of ex- 
cited states corresponding to the first excitations of Na+. 
These excited states lay between the two curves shown in 



3 



•5 15 




Na 2Q + +Na 



10 15 20 25 30 

R(a„) 

FIG. 2: Effective interaction potentials in the harpooning 
model in the case of the of Na^o + Na — > Na^x sticking reac- 
tion. Dotted lines: diabatic potentials. Solid lines: adiabatic 
potentials. The coupling is taken as the hopping integral in 
Na,2 at the distance r. The impact parameter is equal here 
to its critical value b max (see text). 



Figure [2] and might be responsible for quenching towards 
non-reactive channels 
at low collision energy. 



24j . Our model thus applies only 



Practically however, nuclcation 
and evaporation processes involve thermal range collision 
energies which can be considered as small in the sense de- 
fined here. Let us now examine in more details the size 
dependance of a n . First note that there are no strong 
structures. We introduced in our model the experimental 
values for the IP [3(| and for the polarisabilities 31 , 12] • 



Although these values exhibit the well known variations 
due to electronic shells closure [33| . no magic numbers 
clearly appear in the cross sections. This is consistent 
with the unstructured mass spectra observed when clus- 



ters are produced at low temperature [34 , 

The monotonic component does not follow a geomet- 
rical scaling law in n 2 / 3 . Actually, this behavior is not 
really surprising. It has been experimentally observed in 
heterogeneous sticking reactions of small mass selected 
clusters 2l|. The effect on homogeneous cluster growth 
of a discrepancy from the hard sphere approximation had 
already been theoretically suggested in the case of neu- 
tral particles interacting through a Van-Der-Waals inter- 
scales as n 1 / 3 . 



action [17|| . In this scales as '. Surprisingly, 

this kind of effect had never been analyzed nor a fortiori 
experimentally demonstrated in the case of charged clus- 
ters. This non geometrical sticking law has consequences 
for instance on clusters decay rates and kinetic energy re- 
lease distributions (KERD) which are gene rally analyzed 
in the frame of statistical physics models [16[ . Most of 
the models used today are based on the microreversibility 
(or detailed balance) principle, whose basic assumption 
is summarized by the following relation [TBI, [l6| : 
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where 7„ is the cluster decay rate, (3 n -i the sticking rate, 
E n the internal energy, D n the dissociation energy and 
e the kinetic energy of the evaporated atom. fi n , VL atom 
and f2 n _i are the densities of state respectively of the par- 
ent cluster, the evaporated atom and the product cluster. 
(3 n -i is related to the sticking cross section by the rela- 
tion /3„_i = a n -ipv , where p is the volumic density 
of atoms and v the atom velocity. The point we lay the 
stress on here is that in this microreversible approach the 
decay rate j n or the KERD [H, [29| are directly related 
to the sticking cross section <r n _i. 

Dissociation energies are deduced from decay times r„ 
at cluster temperature T through approximations of re- 
lation [3 that can be written to a good accuracy as 



cr„_i A n exp (-D n /(k B T)) 



(3) 



A n is a prefactor independent of a n whose exact expres- 
sion depends on the model. Dissociation energies are de- 
duced from t„ by inverting relation [3] For sodium clus- 
ters, the experimental temperature must be of the order 
of 400K in order to observe dissociation in the typical 
accessible time range |35l [36j . At this temperature, the 
cross section becomes geometrical at about n«200 as de- 
duced from our model (see fig. [I]). One can estimate the 
order of magnitude of the error made by using a geomet- 
rical sticking cross section a geo . Replacing a geo by our 
estimated value at 400K (see fig.Q]) leads to AD = ±4%. 
This error is at least similar to experimental uncertainty 
and should be corrected in future experiments on sodium. 
Unfortunately, a general form of the correction cannot be 
inferred here since it is strongly system-dependent, essen- 
tially through the energy shift and the coupling between 
the two (or more) electronic states involved. 

Nucleation in the frame of the CNT also depends in 
principle on the SCS scaling law. Nucleation theories aim 
to estimate the nucleation rate Ji defined as the number 
of droplets of size i created from gas phase per unit time 
and volume. Although a number of different versions 
have been developed since the work of Becker and Doring, 
the basic assumptions remains the same (a, Iff . CNT is 
derived from the population equations: 



— = ft_mi_i 
ot 



fork + "f l+ in l+1 = Jj_i - Jj(4) 



where and 7 are respectively growth and decay rates. 
The CNT considers the steady state solution such as the 
flux J is constant for all i. This solution Jm is valid 
in a medium time range. Growth at very short delays 
will briefly be examined later. For very large delay, the 
nucleation phenomenon is qualitatively different, notably 
involving coalescence. Moreover, gas depletion prevents 
particles from growing indefinitely. All CNT-based theo- 
ries arrive at an essentially common expression [H, 0, [3^ | : 



Jm = Kfa exp (-AG(**)/(fc B T)) 



(5) 
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where ks is the Boltzmann constant, T the temperature, 
AG(i*) the formation free energy of the critical cluster 
of size i*. The critical size i* is defined by = 0, 

which is equivalent to ^* — 1 [6]. The prefactor K 
does not depend on the cross section. K includes the so- 
called Zeldovich factor [38| which does not depend on the 
sticking scaling law provided that it is a smooth function 
of the size 0, @] • Relation [2l thus the microreversibility 
principle, makes i* independent of the SCS scaling law. 
This had already been pointed out by Vasil'ev and Reiss 
[3^ |. Finally, the nucleation rate in the constant flux 
regime is modified only linearly through . In CNT, 
Pi* is always calculated under the assumption of hard 
sphere collisions. In CUNT, however, many authors con- 
sider that the stronger interaction between charged clus- 
ters and atoms leads to an enhancement of /3j* . They 
introduce a so-called Enhancement Factor (EF) If 
the vapor molecule has no permanent dipolc, EF is gen- 
erally considered as negligible, as far as clusters bear a 
single charge [401 ] - However, when evaluating EF, only 
the interaction between the charged cluster and the neu- 
tral particle is considered. It is demonstrated here that, 
at least in the case of sodium, the harpooning mecha- 
nism strongly modify the SCS's. An example of EF(n) 
at 400 K is readily deduced from figffl it is the ratio 
C400 K/cgeo- EF depends on the system, the tempera- 
ture, the coupling between the electronic states. So we 
cannot derive any general rule for estimating EF. Nev- 
ertheless, we have shown that more attention should be 
paid to this factor. 

Let us focus now briefly on short times. There is a so- 
called lag time th that characterizes the time spent before 
the steady state is reached. The nucleation rate J follows 
approximately the law J(t) ~ Jm(1 — exp (—B0j*t)), 



cluster evaporation measurement results. 
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where the prefactor B do not depend on the SCS [41 1. 
Here the change could be more significant; however, the 
effect on the steady-state nucleation rates is not signifi- 
cant since the lag time b y it self is very short in vapor- 
to-liquid transformations 4l| . Far from equilibrium phe- 
nomena may therefore be more influenced by a change in 
the sticking scaling law. Cluster nucleation in a super- 
sonic expansion is an example of such phenomena. As 
mentioned above, Vigue et al already studied the case 



of neutral clusters [171 ] ■ The final average cluster size 



rif can be represented by empirical scaling laws of the 
type n/ w T 1 ^ for a sticking law in n a . The scaling 
parameter T introduced by Hagena only depends on the 
geometry, the pressure and the temperature in the source 
[421 ] . The example calculated at thermal energy on figure 
[1] shows that a can be close to 0. The change in n/ is 
obviously not negligible. 

In summary, absolute SCS of mass selected sodium 
clusters have been measured for the first time. The non 
geometrical sticking law that depends on the collision 
energy is shown to influence not only nucleation but also 
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